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Abstract — The equivalent circuit model for artificial mag- 
netic materials based on various arrangements of split rings 
is generalized by taking into account losses in the substrate 
or matrix material. It is shown that a modification is needed 
to the known macroscopic permeability function in order to 
correctly describe these materials. Depending on the dominating 
loss mechanism (conductive losses in metal parts or dielectric 
losses in the substrate) the permeability function has different 
forms. The proposed circuit model and permeability function 
are experimentally validated. Furthermore, starting from the 
generalized circuit model we derive an explicit expression for the 
electromagnetic field energy density in artificial magnetic media. 
This expression is valid at low frequencies and in the vicinity 
of the resonance also when dispersion and losses in the material 
are strong. The presently obtained results for the energy density 
are compared with the results obtained using different methods. 

Index Terms — Artificial magnetic materials, permeability func- 
tion, circuit model, energy density 



I. Introduction 

Artificial electromagnetic media with extraordinary proper- 
ties (often called me tamate rials) attract increasing attention 
in the microwave community. One of the widely studied 
subclasses of metamaterials are artificial magnetic materials 
operating in the microwave regime, e.g. [l]-[8]. Broken loops 
have been used as one of the building blocks to implement 
double-negative (DNG) media [9], [10], in addition to this, ar- 
tificial magneto-dielectric substrates are nowadays considered 
as one of the most promising ways to miniaturize microstrip 
antennas [11]-[17]. 

The extraordinary features of metamaterials call for careful 
analysis when studying the fundamental electromagnetic quan- 
tities in these materials. Recently, a lot of research has been 
devoted to the definition of field energy density in DNG media 
[18]-[21]. Authors of [18], [19], [21] derived the energy den- 
sity expression starting from the macroscopic media model and 
writing down the equation of motion for polarization (electric 
charge) or magnetization density in the media. Furthermore, 
complex Poynting theorem was used to search for expressions 
having the mathematically correct form to be identified as 
energy densities. Following the terminology presented in [21] 
we call this method "electrodynamic method" (ED). Tretyakov 
used in [20] another method: Starting from the material mi- 
crostructure, an equivalent circuit representation was derived 
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for the unit cell constituting specific artificial dielectric and 
magnetic media. Lattices of thin wires and arrays of split-ring 
resonators were considered in [20]. The stored reactive energy, 
and, furthermore, the field energy density, were calculated 
using the classical circuit theory. Authors of [21] later called 
this method "equivalent circuit method" (EC). Though these 
two methods apply to media with the same macroscopic 
permeabilities and permittivities, they are fundamentally very 
different, as will be clarified later in more detail. Moreover, 
the final expressions for the field energy density in artificial 
magnetic media given in [20] and [21] differ from each other 
One of the motivations of this work is to clarify the reasons 
for this difference. 

Here we concentrate only on artificial magnetic media 
and set two main goals for our work: 1) To understand the 
differences and assumptions behind the ED and EC-methods 
when deriving the field energy density expressions. We verify 
using a specific example (magnetic material unit cell) that 
in the presence of non-negligible losses one always should 
calculate the stored energy at the microscopic level. 2) To 
generalize the previously reported equivalent circuit repre- 
sentation for artificial magnetic media [20]. The generalized 
circuit model takes into account losses in the matrix material. 
It is shown that this generalization forces a modification to the 
widely accepted permeability function used to macroscopically 
describe artificial magnetic media. The generalization has 
a significant importance as it is shown that in a practical 
situation matrix losses strongly dominate over conductive 
losses. The proposed circuit representation and permeability 
function are experimentally validated: We measure the mag- 
netic polarizability of a small magnetic material sample and 
compare the results with those given by the proposed analytical 
model and the previously used model. The results given by 
the proposed model agree very well with the measured results, 
whereas the old model leads to dramatic overestimation for the 
polarizability. Using the generalized circuit model we derive 
an expression for the field energy density in artificial magnetic 
media. This expression is compared with the results obtained 
using the ED-method in [21]. Reasons for the differences in 
the results are discussed. 

II. Electrodynamic method vs. Equivalent circuit 

METHOD 

It is well explained in reference books (e.g. [22], [23]) that 
for the definition of field energy density in a material having 
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non-negligible losses one always needs to know the material 
microstructure. First of all, the reactive energy stored in any 
material sample is a quantity that can be measured. When 
the material is lossless, no information is needed about the 
material microstructure for this measurement. Indeed, we can 
measure the total power flux through the surface of sample 
volume and, since there is no power loss inside, we can use the 
Poynting theorem to determine the change in the stored energy. 
This is the reason why the field energy density in a dispersive 
media with negligible losses can be expressed through the 
frequency derivative of macroscopic material parameters. In 
the circuit theory, the same conclusion is true for circuits that 
contain only reactive elements: It is possible to find the stored 
reactive energy in the whole circuit knowing only the input 
impedance of a two-pole [24]. 

Simple reasoning reveals, that in the presence of non- 
negligible material losses the above described "black box" 
representation and direct measurement are not applicable: 
Without knowing the material microstructure or the circuit 
topology we do not know which portion of the input power 
is dissipated and which is stored in the reactive elements. 
Thus, the energy stored in a lossy media cannot be uniquely 
defined by only utilizing the knowledge about the macroscopic 
behavior of the media [22], [23]. 

When defining the energy density in a certain material using 
the ED-method one first writes down the equation of motion 
for charge density or for magnetization in the medium using 
the macroscopic media model [19], [21]. We stress that this 
equation is the macroscopic equation of motion, containing 
the same physical information as the macroscopic permittivity 
and permeability. Further, complex Poynting theorem is used 
to identify the mathematical form of the general macroscopic 
energy density expressions. Having the form of these expres- 
sions in mind one searches for similar expressions in the 
equation of motion and defines them as energy densities. 
The problem of the ED-method method is the fact that it 
only utilizes the knowledge about the macroscopic behavior 
of the media, which, as explained above, in not enough. 
The aforementioned difficulty is avoided in the EC-method 
[20]. Based on the known microscopic medium structure, 
one constructs the equivalent circuit for the unit cell of the 
medium. Careful analysis is needed to make sure that the 
circuit physically corresponds to the analyzed unit cell. After 
this check the stored reactive energy and the corresponding 
field energy density can be uniquely calculated using the 



Fig. 1. Split-ring resonator loaded with an infinitesimally small lumped 
circuit. 



classical circuit theory. 

Next, we illustrate the difference between the ED- and 
EC-methods using a specific case of split-ring composites. 
Consider the split-ring resonator (SRR) shown in Fig. ^ 
Following a example given in [23], we load the SRR with an 
electrically infinitesimally small circuit consisting of lumped 
elements. Let us assume that the additional inductance and 
capacitance are chosen to have values = tR^,C^ = 
t/R^,t > 0. A simple check reveals that in this case the 
input impedance of the load circuit is frequency independent 
and purely resistive: Zi„ ~ R^. When the SRR is electrically 
small, it can be represented as a resonant contour and the 
total loss resistance reads i?tot = R + R^, where R is the loss 
resistance due to the finite conductivity of the loop materials. 
Let us further assume that the ring is made of silver and the 
value of R^ is chosen so that i?tot is the same as R for an 
unloaded SRR made of copper. In this case the macroscopic 
descriptions of a medium formed by the loaded SRRs made 
of silver with additional loads and unloaded SRRs made of 
copper are exactly the same. Thus, the ED-method predicts 
the same value for the reactive energy stored in these two 
media. Inspection of Fig. ^ clearly indicates, however, that 
this is not the case. There is additional energy stored in the 
load inductance and capacitance, which is invisible on the level 
of the macroscopic permeability description. Proper definition 
of the stored energy must be done at the microscopic level, 
which is possible with the equivalent circuit method. 

III. Equivalent circuit method: Brief revision of 

EARLIER RESULTS 

A commonly accepted permeability model as an effective 
medium description of dense (in terms of the wavelength) 
arrays of split-ring resonators and other similar structures reads 



(1) 



oJq — uj^ + jajT 

(see e.g. [1], [2], [4], [8].) Above, A is the amplitude factor 
(0 < A < 1), llIq is the undamped angular frequency of the 
zeroth pole pair (the resonant frequency of the array), and F is 
the loss factor The model is obviously applicable only in the 
quasi-static regime since in the limit w — > oo the permeability 
does not tend to /iQ. At extremely high frequencies materials 
cannot be polarized due to inertia of electrons, thus, a phys- 
ically sound high frequency limit is hq [22]. However, {TJ 
gives correct results at low frequencies and in the vicinity of 
the resonance. This is the typical frequency range of interest 
e.g. when utilizing artificial magneto-dielectric substrates in 
antenna miniaturization [13], [16], [17]. The other relevant 
restriction on the permeability function is the inequality [22] 



dijj 



(2) 



valid in the frequency regions with negligible losses. Physi- 
cally the last restriction means that the stored energy density 
in a passive linear lossless medium must be always larger than 
the energy density of the same field in vacuum. Macroscopic 
model ([0 violates restriction 1 at high frequencies, which 
is another manifestation of the quasi-static nature of the 
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Fig. 2. a) Magnetic material sample in the probe magnetic field of a tightly 
wounded long solenoid, b) Equivalent circuit model, losses in the matrix 
material are not taken into account. 
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Fig. 3. a) Magnetic material sample in the probe magnetic field of a tightly 
wounded long solenoid, b) Equivalent circuit model taking into account losses 
in the matrix material. 



model. In the vicinity of the magnetic resonance the effective 
permittivity of a dense array of spht-ring resonators is weakly 
dispersive, and can be assumed to be constant. 

In [20] the energy density in dispersive and lossy magnetic 
materials was introduced via a thought experiment: A small 
(in terms of the wavelength or the decay length in the 
material) sample of a magnetic material [described by Q] was 
positioned in the magnetic field created by a tightly wounded 
long solenoid having inductance Lq, Fig. |2^. The insertion 
changes the impedance of the solenoid to 



Z{uj) = jujLQfi^{uj) = jujLo + 



UJ, 







(3) 



'2 J- jujY 

The equivalent circuit with the same impedance was found to 
be that shown in Fig. |2]3 [20] with the impedance seen by the 
source 



Z{uj) = jiuLo + — 



LC 



which is the same as (|3j if 



LLo 



1 

LC 



^0' 



L 



(4) 



(5) 



The aforementioned equivalent circuit model is correct from 
the microscopic point of view since the modeled material is a 
collection of capacitively loaded loops magnetically coupled 
to the incident magnetic field. An important assumption in 
[20] and in the present paper is that the current distribution 
is nearly uniform over the loop. This means that the electric 
dipole moment created by the exciting field is negligible as 
compared to the magnetic moment. The electromagnetic field 
energy density in the material was found to be [20] 



(wg _ ^2)2 ^ ^2r2 



\H\ 



(6) 



In [20] only losses due to nonideally conducting metal of 
loops were taken into account, and losses in the matrix mate- 
rial (substrate material on which metal loops are printed) were 
neglected. It will be shown below that neglecting the matrix 
losses can lead to severe overestimation of the achievable 
permeability values. 

IV. Generalized equivalent circuit model and 

PERMEABILITY FUNCTION 

Losses in the matrix material (typically a lossy dielectric 
laminate) can be modeled by an additional resistor in parallel 



with the capacitor. Indeed, if a capacitor is filled with a lossy 
dielectric material, the admittance reads 



Y = jujC{e' - je") ^ jujCe' + ujCe" 



(7) 



where the latter expression denotes a loss conductance. Thus, 
the microscopically correct equivalent circuit model is that 
shown in Fig. |3j>. The impedance seen by the source can be 
readily solved: 

. , , jc.H-P/L + u;^MyiLCR^) 

Zj — ILUL/n -\ 7z — . (Q) 

The macroscopic permeability function corresponding to this 
model reads 

V (l + l!7)ra-"'+J"(T+c7i7)/ 
Comparing Q and (|5Jl we immediately notice that Q is an 
insufficient macroscopic model for the substrate if the losses in 
the host matrix are not negligible. A proper macroscopic model 
correctly representing the composite from the microscopic 
point of view is 

/ Au'^ - jujB \ 

Equation ^ is the same as dlOt if 



M2 



A 



Af2 



LLq LLqC Rd 



CRd 



R \ 2 —2 



(11) 



Above we have denoted o^q = 1/{LC). The macroscopic 
permeability function of different artificial magnetic materials 
can be conveniently estimated using dlOt . as several results 
are known in the literature for the effective circuit parameter 
values for different unit cells, e.g. [2], [6], [8]. 

For the use of dlOl l it is important to know the physical 
nature of the equivalent loss resistor i?d- If losses in the matrix 
material are due to finite conductivity of the dielectric material, 
the complex permittivity reads 



( = ( - je = e 



(12) 



where a is the conductivity of the matrix material. Thus, we 
see from that the loss resistor is independent from the 
frequency and can be interpreted as a "true" resistor Moreover, 
in this case the permeability function is that given by dlOt . 
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However, depending on the nature of the dielectric material the 
loss mechanism can be very different from (I12> . and in other 
situations the macroscopic permeability function needs other 
modifications. For example, let us assume that the permittivity 
obeys the Lorentzian type dispersion law 

e = e'(l + — S ^1, (13) 



where lo'^ is the angular frequency of the electric resonance, 
C is the amplitude factor and A is the loss factor Moreover, 
we assume that the material is utilized well below the electric 
resonance, thus, lo ^ lo'^. With this assumption the permittivity 
becomes 

('{1 + C)- jLoe'CK/io'^. (14) 

We notice from (0 that in this case the equivalent loss resistor 
i?d becomes frequency dependent: 



i?d oc 



(15) 



and the permeability function takes the form 



(16) 

where K is a real-valued coefficient depending on the dielec- 
tric material. For other dispersion characteristics of the matrix 
material the permeability function can have other forms. 

V. Experimental validation of the proposed 

CIRCUIT model and PERMEABILITY FUNCTION 

In this section we present experimental results that validate 
the generalized equivalent circuit model and the corresponding 
macroscopic permeability function. The measurement cam- 
paign and the experimental results are described in detail in 
[8]. For the convenience of the reader we briefly revise the 
main steps of the measurement procedure. 

The measured artificial magnetic particle, metasolenoid, 
is schematically presented in Fig. 0^. In [8] the effective 
permeability of a medium densely filled with infinitely long 
metasolenoids was derived in the form 



Ztntd 



(17) 




a) 



b) 



Fig. 4. a) A schematic illustration of the metasolenoid. b) A finite-size 
metasolenoid approximated as a magnetic ellipsoid. 



where V^. is the volume filling ratio, S = a x b is the cross 
section are of the ring, d is the distance between the rings, 
and the total effective impedance was presented in the form 



1 



off- 



(18) 



off 



For the experimental validation a finite-size metasolenoid 
was approximated as an ellipsoid cut off from a magnetic 
media described by ( I17> . Fig. The magnetic polarizability 
of the ellipsoid amm was analytically calculated using the 
classical mixing theory [25]. The scattering amplitude of 
an electrically small material sample was measured using a 
standard parallel plate waveguide, and the theory of waveguide 
excitation was used to calculate the magnetic polarizability 
from the measured results. It is worth noting that the magnetic 
polarizability is a function of the magnetic permeability due to 
the mixing process. On the other hand, permeability is defined 
using the equivalent circuit. Fig. |3j3. Thus, the magnetic 
polarizablity of the measured sample contains all the relevant 
data for validating both the proposed equivalent circuit and 
the permeability function. 

Though it is not explicitly mentioned in [8], substrate losses 
were taken into account when analytically calculating the 
magnetic polarizability of the sample. The authors used ( I18t 
to define the total impedance of the metasolenoid unit cell, 
however, complex permittivity was used when calculating the 
effective capacitance. Thus, the equivalent circuit of the unit 
cell used to analyze the measured sample is the proposed 
circuit shown in Fig.|3j5. It can easily be verified using the data 
presented in [8] that the following expression for the effective 
impedance (derived using the circuit in Fig.lSj') exactly repeats 
the analytical estimation for the magnetic polarizability of the 
sample: 



^tot 



jwLcff + R, 



Rd 



•off 



1 + juiCcsRd ' 



(19) 



where Rd — e' / [e" LuCcff)- The analytically calculated [2'tot 
given by M9\ is used in ( I17H and measured magnetic polar- 
izabilities are shown in Fig. [S] The measured and calculated 
key parameters are gathered in Table |l] The polarizability and 
permeability values in Table U are the maximum values. 

The measured results agree rather closely with the analytical 
calculations when the proposed model is used. The slight 
difference in the resonant frequencies, and the slightly lowered 
polarizability values in the measurement case are most likely 
caused by limitations in the accuracy of the manufacturing 
process: The implemented separation between the rings is 
probably slightly larger than the design value. This lowers 
the effective capacitance and is seen as a weakened magnetic 
response and the higher resonant frequency. Moreover, the 
measured Im{ainm} clearly indicates that the effect of the 
lossy glue used to stack the rings is underestimated in the 
analytical calculations. In the analytical calculations tan (5 — 
0.002 was used for the total loss tangent [8]. A loss tangent 
value tan (5 = 0.0025 would accurately produce the measured 
polarizability values, and the bandwidths (defined from the 
Imlamm} curve) would visually coincide. 
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TABLE I 

Analytically calculated and measured properties of the 
metasolenoid sample. 
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1.09 
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4000 
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If the matrix losses are neglected in the circuit model 
[i?d — > oo in (I19H. the analytical calculations lead to dra- 
matically overestimated polarizability and permeability values. 
It is therefore evident that the proposed generalization of the 
circuit model and the permeability function have a significant 
practical importance. Though the model has been validated 
using a specific example, we can conclude that matrix losses 
can strongly dominate over conductive losses in structures 
where the unit cells are closely spaced. This is physically well 
understandable since in this case most of the flux is forced 
inside the substrate. 

VI. Electromagnetic field energy density 

Following the approach introduced in [20] we will next 
generalize the expression for the energy density in artificial 
magnetics using the experimentally validated circuit model. 
In the time-harmonic regime the total stored energy reads 
(notations are clear from Fig. (Sji) 
1 



W^-{Lo\I\^ 



C\VcV) 



[(i + f)xb 



(20) 



(21) 
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Fig. 5. The analytically calculated (proposed model) and measured magnetic 
polarizability. 



Using the notations in M 11 1 the stored energy can be written 
as 



w^\l,\i\^ 



1 + 



(22) 



The inductance per unit length of a tightly wound long 
solenoid is Lq = l^on^S, where n is the number of turns 
per unit length and S is the cross section area. The relation 
between the current / and magnetic field H inside the solenoid 
is / = H/n. Thus, the stored energy in one unit-length section 
of the solenoid reads 

2 



1- 



w2)2+c^2(r + rd)2 

(23) 

from which we identify the expression for the electromagnetic 
field energy density in the artificial material sample: 



f~2 



\H\ 



(24) 



We immediately note that if there is no loss in the matrix 



material (i?d - 
reduces to (|6j. 



oo and Va 



0), then Lo^ 



and ( l24l i 



A. Comparison with the results obtained using the ED-method 

The above derived result differs from the result found in 
[21] using the ED-method: 



t^2[(,^2 _^2)2 +^2r2] 



H 



(25) 



The procedure and the underlying assumptions to obtain ( I25t 
have been briefly revised in Section II. The classical expression 
for the magnetic field energy density in a media where 
absorption due to losses can be neglected reads [22], [23] 



H 



(26) 



2 duj 

It is seen that in the presence of negligible losses [F — > in 
Q] the energy density result given by M5\ is the same as the 
result predicted by the classical expression \2bl . However, MAI 
predicts a different result. Authors of [21] use this fact to state 
that the result obtained using the ED-method is more internally 
consistent than the result obtained using the EC-method. 

The EC-method is known to give a perfectly internally 
consistent result for the energy density in dielectrics obeying 
the general Lorentzian type dispersion law [20]. The general 
Lorentz model is a strictly causal model. This is, however, 
not the case with the modified Lorentz model Q. As is 
speculated already in [20], the reason for the difference in 
the results obtained using ( I24> and M6\ in the small-loss limit 
is related to the physical limitations of the quasi-static perme- 
ability model Q. Thus, when ([ij is used as the macroscopic 
media description, MA\ should be used also in the presence 
of vanishingly small losses. The ED-method, though being 
internally consistent with the classical expression, predicts 
unphysical behavior at high frequencies: At high frequencies 
the energy density given by the ED-method is smaller than 
the energy density in vacuum (when there are no losses, this 
unphysical behavior takes place at frequencies u > ^/i^JJo, 
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Fig. 6. Magnetic field energy density given by different expressions, ojq is 
the resonant frequency of tlie metasolenoid medium. 

where restriction 1 is violated). Fig. la and lb in [21]. This 
behavior is avoided with the result obtained using the EC- 
method, since that approach is based on the microscopic 
description of the medium, which is always in harmony with 
the causality principle. 

Fig- m depicts the normalized magnetic field energy density 
in a medium formed by the metasolenoids introduced in the 
previous section. The amplitude factor A = 1, and the loss 
factors have been estimated using (lll> and the data presented 
in [8]. In this particular example the energy densities given by 
i24\ and (|6jl give practically the same result over the whole 
studied frequency range (the result given by (|6j is not plotted 
in Fig. |6l since the result visually coincides with the EC- 
result). This is due to the fact that large values of ut and 
ojo mask the effect of F and Fd in i24i and The results 
given by the ED-method and the classical expression i26\ 
also visually coincide. However, as was mentioned above, the 
energy density expression given by the ED-method predicts 
the same nonphysical behavior as the classical expression: 
The field energy density is smaller than the energy density 
in vacuum at frequencies uj > a/Swq- 

Conclusion 

In this paper we have explained differences between recent 
approaches used to derive field energy density expressions 
for artificial magnetic media. The equivalent circuit model 
of split-ring resonators and other similar structures has been 
generalized to take account losses in the dielectric matrix 
material. It has been shown that a modification is needed to the 
macroscopic permeability function commonly used to model 
these materials in the quasi-static regime. Moreover, depending 
on the nature of the dominating loss mechanism in the matrix 
material, the permeability function has different forms. The 
proposed circuit model and the modified permeability function 
have been experimentally validated, and it has been shown that 
in a practical situation matrix losses can dramatically dominate 



over conductive losses. Using the validated circuit model we 
have derived an expression for the electromagnetic field energy 
density in artificial magnetic media. This expression is valid 
also when losses in the material cannot be neglected and 
when the medium is strongly dispersive. The results have been 
compared to the recently reported results. 
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